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Chapter 1 
Introduction 


Rendezvous of chaser and target vehicles has been and will continue to be an important 
technological capability required to support space exploration and space utilization 
programs. This ability to safely and effectively meet up with another object in space is 
necessary for various types of missions, from transfer of personnel or supplies to rescue, 
retrieval, repair, inspection, or assembly of orbital units. Rendezvous capability made the 
Apollo moon landings possible and allows the Space Shuttle to retrieve and repair Earth 
satellites. It will be required for the build up and resupply of Space Station Freedom and 
will be essential for the development of Lunar and Martian outposts. 


The mission planner must select the best rendezvous scheme under diverse, and possibly 
conflicting, mission constraints. For example, the planner must trade off fuel expenditure 
against duration of the rendezvous phase. Throughout rendezvous, the relative range 
from one craft to the other should be small enough to allow for navigation. The intervals 
between the various burns must be large enough to permit tracking and to allow 
preparation for the subsequent maneuver. For safety and control purposes, the closing 
velocity should not be excessive. The final approach geometry should produce lighting 
conditions which are favorable to sensor tracking. Furthermore, certain constraints may 
be imposed on the maneuvers to allow for manual modes of operation if the automatic 
systems fail. 
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1.1 Background 


Rendezvous technologies developed for the Apollo and Skylab programs involved chaser 
and target vehicles in near-circular orbits (orbits of very low eccentricity). Both 
programs employed rendezvous profiles which consisted of similar maneuvers to control 
the closure phasing and height of the active vehicle (chaser) relative to the passive vehicle 

(target). 

For typical Apollo and Skylab missions (Figure 1.1), the first few maneuvers placed the 
chaser in a coplanar intermediate phasing, or staging, orbit. This orbit was coelliptici 
with that of the target at a slightly lower altitude so that the chaser caught up slowly from 
below and behind. These coelliptic orbits provided minimal altitude variation as the 

chaser gradually closed on the target. 



Figure 1.1 Apollo/Skylab Coelliptic Rendezvous Profile 


iThe product of the semi-major axis and the eccentricity of the staging orbit is equal to that of the target 
orbit, and the line of apsides of the staging orbit is coincident with that of the target orbit. 
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Once the chaser had reached a desired elevation angle of the line-of-sight from the chaser 
to the target (in the chaser's local vertical frame), the chaser would execute the intercept 
maneuver and place itself on an intercept trajectory (the transfer phase). Midcourse 
correction maneuvers were made to maintain the intercept trajectory. The final braking 
maneuvers were then executed to provide the proper attitude and rate of closure for 
docking. 

This rendezvous strategy allowed for standardization of transfer profiles. Use of the 
coelliptic phasing orbits resulted in reasonably low, nearly constant closing rates prior to 
transfer phase initiation (TPI) and standard lighting and approach conditions at 
completion. This simplified training and mission planning and allowed backup 
calculation of the burn, independent of the onboard computer. In addition, small bum 
magnitudes were required, minimizing the effects of a poorly executed maneuver. While 
the use of phasing orbits requires more fuel and takes longer than a direct transfer, the 
benefits far outweigh the costs. 

The Shuttle rendezvous procedure is somewhat different in that the phasing portion of the 
profile does not maintain a coelliptic closure with the target, but instead employs an orbit 
which intercepts the target altitude (Figure 1.2) allowing for "stable" active vehicle 
standoff conditions if the need arises. Transfer is initiated from this staging orbit. After 
transfer is initiated (TI), the second correction maneuver (MC2) employs the same 
elevation angle condition as used in Apollo and Skylab. These maneuvers result in 
closure rates and relative geometries similar to those in Apollo and Skylab rendezvous. 
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Figure 1.2 Shuttle Rendezvous Profile 
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While the Apollo and Shuttle rendezvous strategies were developed for targets in near- 
circular orbits, missions to Mars may require rendezvous with targets in non-circular 
orbits. Some of the current mission scenarios have created an interest in the development 
of highly elliptic rendezvous capability. In these scenarios, the heavy interplanetary 
vehicle is parked in a high-energy, highly elliptic orbit instead of the traditional low- 
energy circular orbit. Only the small landing craft is taken down to a low circular orbit 
and/or down to the surface. This approach avoids the fuel penalties associated with 
taking such a heavy platform down to a low altitude orbit only to have to lift it back up 
later. 

1.2 Error Sources 

Most rendezvous studies have focused on fuel minimization with little regard for 
operational or navigation concerns. With the fuel saved by leaving the heavier platform 
in orbit, rendezvous studies can now focus on other criteria. Obviously, the result must 
still be reasonably fuel efficient and cannot be done with a total disregard to fuel usage. 

In general, applying an incorrect bum at transfer initiation - whether because of targeting 
errors due to navigation uncertainties or because of maneuver execution errors - will 
result in a failed rendezvous. Since the primary goal for rendezvous missions is to bring 
two vehicles to intercept, the actual time that intercept occurs is of less importance than 
the fact that it does occur at some point in the trajectory. 

Any error in bum estimation or implementation will result in an error in the position of 
the chaser at the nominal intercept time. Angular measurements supply direct information 
about line-of-sight (LOS) from the chaser to the target, but only very indirect information 
about range along that LOS. On the other hand, with range only measurements, the 
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chaser will have very good knowledge of its range to the target and imperfect knowledge 
of the direction of the LOS. Furthermore, the beginning and ending of a bum can only be 
controlled to the frequency at which the onboard software is cycled. Also, translation 
engines have a finite rise and tail off time, further affecting the Av imparted to the chaser. 
In addition, control of the burn direction depends on the accuracy of the inertial 
measurement unit (IMU), which drifts over time, as well as the attitude deadband of the 
active vehicle. By allowing a small variation in the transfer time of the intercept 
trajectory, we gain a much-needed degree of freedom that can be used to compensate for 
selected one-dimensional errors in estimation or execution of the intercept maneuver (TPI 
in Figure 1.1). 

1.3 Overview 

This thesis will develop rendezvous maneuver point conditions which provide error 
tolerance. Targeting of an intercept maneuver is only as accurate as estimate of the 
relative state between the chaser and the target. Even if the relative state were known 
exactly (so that the required change in velocity can be calculated exactly), intercept bum 
execution errors would still limit the ability to rendezvous successfully. 

1.3.1 Time as a Degree of Freedom 

Given that it is more important for rendezvous to take place at some point along the 
trajectory than that it occur at a particular time, we can allow the time of intercept to vary 
in order to ensure that a rendezvous does, in fact, occur. There are two ways to view the 
use of this degree of freedom to yield insensitivity to errors. We can view the effect of 
the one-dimensional errors as a velocity error at transfer initiation and search for times 
when this velocity error is aligned with the one direction of freedom in post-bum 
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velocity. Or we can see the effect of the errors as position errors at intercept and perform 
the same type of alignment with the relative velocity of the two vehicles at intercept. 
This thesis will take the first approach. The second approach results in equivalent 
maneuver point conditions, but in terms of the Kepler transition matrix instead of those of 
the Lambert transition matrix. These alternate forms are developed in Appendix A. 
Chapter 2 identifies the freedom in the required burn obtained by allowing small 
variations in time of intercept. 

1.3.2 Execution Error Tolerance 


Chapter 3 presents the analytic development of maneuver point conditions which will 
result in a rendezvous which is insensitive to errors in maneuver execution (i.e., bum 
magnitude, direction, and ignition time). We will demonstrate the validity of these 
maneuver point conditions by applying them to rendezvous with targets in both circular 
and elliptic orbits. We will first place the target in a circular Earth orbit at 350 km with 
the chaser in a cocircular orbit 20 km below, chosen to be representative of Apollo and 
Skylab missions. 2 Because of the interest in highly elliptic orbits for missions to Mars, 
we will then place the target in a highly elliptic Martian orbit with a period 
approximately equal to one Martian day (r p = 500 km and r a = 33000 km) and the chaser 
in either a scaled elliptic^ or coelliptic orbit, (Figure 1 .3) each with 20 km decrease in the 
semi-major axis. 


2 These robust maneuver point conditions arc also readily applicable to transfer from the more complex 
staging orbits used by the Shutde. 

•*The eccentricity of the staging orbit is equal to that of the target orbit, and the line of apsides of the staging 
orbit is coincident with that of the target orbit. 
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Figure 1.3 Scaled Elliptic and Coelliptic Staging Orbits 

1.3.3 Navigation Error Tolerance 

In Chapter 4, we will investigate errors in the estimate of the relative state of the chaser. 
We will then use the freedom in intercept time to develop maneuver point conditions 
robust to these types of errors. We will perform deterministic flyout from a maneuver 
point and plot the resulting relative motion as we did in Chapter 3. We will further 
demonstrate the validity of the maneuver point conditions by generating dispersed states 
typical of line-of-sight and range measurements respectively, propagating the states to 
intercept, and observing the focusing effect. 4 


4 The dispersed chaser states converge on the target, although at different times. 



Chapter 2 

Variation in Time of Intercept 


2.1 Derivation of Degree of Freedom 

By allowing a small variation in the time of intercept ( t [), we create a degree of freedom 
that can be used to compensate for any one-dimensional error. Such errors may result 
from bum execution errors, navigational uncertainties, or some combination of both. In 
order to make use of this freedom in initial burn, we must first identify it. The velocity 
that the chaser must achieve at transfer initiation is a function of the transfer time 
involved ( At = t^ -r 0 ), the initial position of the chaser (r c (r 0 )), and final position of the 
target (r T (t x )). 
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This is essentially Lambert's problem: 


v^ = y^(Ar,r c (r 0 ),r 7 .(r 1 )) 


( 2 . 1 ) 


The required bum is equal to the difference between the required velocity and the pre- 
bum velocity of the chaser: 


Av rt<t — y r , <J (Ar,r c (r 0 ),r J -(r 1 )) y c (r 0 ) 


( 2 . 2 ) 


If we allow small perturbations in any of these variables, the burn required for 
rendezvous will change. Taking the full differential of Equation 2.1, 


5(Av^) = %!-5A/+ -pa- Sr c (t 0 )+ «&<«,)-*£«,) 
’ dAt dr c (t 0 ) dr r (t ,) 


(2.3a) 


Or more compactly, 


<5( Av rtq ) = X 0 5 At + Li Sr c {t 0 ) + U Sr T (r, ) - 8)T C (r 0 ) 


(2.3b) 


where from Lambert's problem 


A. 


_ dv rtq 

dAt 


and Lj and L 7 are components of the Lambert "transition" (or sensitivity) matrix 


'Sv c (t Q y 

=o, 

dr c (t 0 ) 

_<5y c (r,) 

L 

_Sr T (t l ) 


■L, L,' 

'5r c (r 0 )' 

l 3 l 4 

_5r 7 (r, ) 


(2.4) 
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r _ dv c (l 0 ) 

dr c (t o) 

r _ §Yc<A) 
dr c (t o) 


r - ^Vc(^o) 
<?/>(/,) 
L _ dv c (r,) 
4 *&(*) 


Now, if we allow the time of intercept to slip by a small amount, SAt = St , . The initial 
position and velocity of the active vehicle do not change, Sr c (t 0 ) = Svc(t 0 ) = 0, but final 


position of the passive vehicle does: 


<5r r ('i) = V7.('i)<5'i 


(2.5) 


The change in the required burn caused by this small time slip is then 


S(Av rtq ) = [A 0 + L 2 y r (r, )]<5r. 


(2.6a) 


The first term on the right hand side of this equation is the change in required bum due to 
a change in intercept time for rendezvous with a fixed point in space. The second term is 
the change in required burn caused by the movement of the target. Together these effects 
determine the direction in which the initial post-burn velocity is free to vary. 

For convenience, we will define 
vv,, = A 0 + L 2 y r (r 1 ) 

so that 

^(Av„,) = w, i 5r 1 (2.6b) 
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Thus, small bum errors which are parallel to w, t result in small final position errors which 
are parallel to the relative velocity vector at rendezvous, and these errors are absorbed by 
small variations in transfer time. Small errors in initial post-bum velocity which are 
perpendicular to vv (i result in small final position errors which are perpendicular to the 
relative velocity vector and therefore result in pure miss distance. 


2.2 Linearity Restrictions 

The derivation includes assumptions of linearity which place restrictions on the 
application of this freedom. First, we assume that the variation in the transfer time is 
small. 


Ar 

and that the change in the position of the target is small 

1^-r (fi )[ j 

Mol 

We also assume that the motion of the target at the intercept point is linear. This 
assumption of linear motion at intercept is often the most restrictive. The length of time 
for which this is true may be quite small, particularly near periapse. When we use this 
degree of freedom to absorb errors, we will have to monitor changes in intercept time to 
ensure that all of these assumptions hold true. 
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Chapter 3 

Execution Error Tolerance 


Any error in implementation of the burn will result in an error in the position of the 
chaser at the nominal intercept time. We will examine three types of bum error: 

magnitude error, in-plane directional error, and ignition time error. The bum can be 
expressed as Av c (r 0 ) = Vc(r 0 )-yc(f 0 ). When the burn is executed exactly, v£(r 0 ) = v rtq 
and Ay c (r 0 ) = Ay re<r 


As Figure 3.1 shows, an error in burn magnitude is parallel to Ay„ fl , while a small error 
in bum direction is perpendicular to Ay rei? . These errors are the two in-plane components 

of executing the wrong burn at the right time. The conditions that result in a rendezvous 
which is robust to small errors in magnitude or direction are developed together in 
Section 3.1. The condition for robustness to an ignition time error (executing the right 
bum at the wrong time) is developed in Section 3.2. 


Bum Magnitude Error 


5(Ai c (f 0 )) l( = £(Ai r< ,) 


A^ 


1 ^ Burn Direction Error 

I <5( A jfc- (r 0 )) j. = ( Ax rtq x i *) 50 


Figure 3.1 Burn Magnitude and Direction Errors 
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3.1 Magnitude and Direction Errors 


Consider the case in which the burn is applied at the correct time but with an error in 
either the magnitude or the direction. Since there is no deviation in initial position or 
velocity of the chaser (Sr c (t 0 ) = Sv c (t 0 ) = 0) and the final position of the target is 

unchanged, the required bum is unchanged. 

In order for rendezvous to occur, the error in the burn must be equal to the acceptable 
deviation in the bum associated with the variation in intercept time: 

<5(Av) = 5(Av r , q ) = vy (i St, (3- 1 ) 

In the case of a burn magnitude error, <5(Ay) l|; 

S(Ay)„ = eAy r<<f = w, t St, ( 3 - 2 ) 


For this to be true, the two vectors must be parallel or anti-parallel: 
Ay„, x vy (i = 0 


(3.3) 


Thus, Equation 3.3 is the maneuver point condition defining insensitivity to bum 
magnitude errors. The corresponding intercept time variation is 


St, 


= e A-««_=!l _ ±£ 




(3.4) 


In the case of an error in burn direction, <5(Ay) x , 
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<5(Av) ± = (Av r< , x i h )S6 = vv (| St, 


(3.5) 


where i h is a unit angular momentum vector (perpendicular to the orbital plane). For this 
to occur, the bum direction error must again be parallel or anti-parallel to the acceptable 
deviation in the bum 


<5(Av) x x vv (i =0 

Av„ ? vv,, =0 (3.6) 


Equation 3.6 is the robust maneuver point condition for burn direction errors. The 
resulting variation in the time of intercept is 




(A /„)•*,,, 

w. • w. 

— *1 — ‘l 


Av 

se = ±-^sd 

w, 


(3.7) 


These two maneuver point conditions are complementary and span the space of in-plane 
bum execution errors. Clearly, these conditions cannot occur simultaneously. A 
maneuver point condition can be chosen which is insensitive to a bum magnitude error or 
one which is insensitive to a bum direction error, but not one which is insensitive to both 
at the same time. However, if a specific combination of magnitude and direction error 
(an error confined to one dimension) is expected, we can develop a maneuver point 
condition which is robust to that error instead. 
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3.2 Ignition Time Error 


If we ddsy execution of the bum by some small amount of time, St 0 , the chasers initial 
state, both position and pre-bum velocity, will change slightly: 


5r c (f o) — HcOo)^o 

Sv~ c (t 0 ) = g (t 0 )St 0 = — p C c (t 0 )St 0 

LM\ 


(3.8) 


Because of this delay, the chaser has less time to arrive at the nominal intercept point: 


SAt = -<5r 0 < 3 - 9 > 

Substituting Equations 3.8 and 3.9 into Equation 2.3, the deviation in the required bum is 

<5(Av„, ) = [- A o + Li Yc (f 0 ) - g c (f 0 )] &o (3 ‘ 10a) 

The first term on the right hand side of this equation is the change in required bum with a 
change in time of execution for rendezvous from a fixed state in space. The next two 
terms are the change in required burn caused by the change in the initial state of the 
chaser. Together these effects determine the direction in which the initial required bum 
varies. 

For convenience, define 

vv, 0 = - A o + Li v'c (t 0 ) - g c (r 0 ) 
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so that 


S(Av rtq ) = w, o 8t 0 (3.10b) 


In order for rendezvous to occur, this deviation in required bum combined with the 
deviation in required bum associated with a deviation in time of intercept must be zero. 

S(Av req ) = vv (i <5r, + vv, o 5r 0 =0 (3.11) 


Equation 3.12 is the maneuver point condition which provides insensitivity to ignition 
time errors 


w, x w, =0 

— * o — *i — 


The corresponding variation in time of intercept is 


w, ■ w, „ 

&, = - ■ 1 = ± 


w , • w. 

'l M 


w, 


w r 


Sl 


(3.12) 


(3.13) 


A delay (executing the right burn at the wrong time) results in a burn error which is 
generally some combination of error in magnitude and direction. In special cases, this 
maneuver point condition may occur simultaneously with one or the other of the 
previously developed conditions. In other words, we may be able to choose a maneuver 
point which is robust to both burn magnitude errors and burn ignition time errors or one 
which tolerates both burn direction errors and burn ignition time errors. 
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3.3 Application of Burn Error Tolerance Conditions to Circular Orbits 


With knowledge of the target's orbit and of the chaser's orbit, the various possible transfer 
alternatives can be searched for those which result in a maneuver point which meets the 
desired condition. If such maneuver points exist, they can be used to achieve tolerance to 
whichever type of error is of most concern. 

For a given circular orbit pair, a nominal rendezvous is completely defined by the relative 
true anomaly of the target ( fyJJo) = 0 ) - Mo)) and die central transfer angle (0). 



For each choice of 0, we will vary f r , (r 0 ) 1 and plot the angle <p (Figure 3.3) between 

the required bum and the change in the required burn caused by the variation in intercept 
time (Figure 3.4). 


1 % (?o) 


must be small to keep fuel consumption reasonable. 
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Figure 3.3 Angle Between Required Burn and Change in Required Bum Insensitivity 
Direction 



When <p = nn, the condition for burn magnitude error tolerance (Equation 3.3) is 
satisfied. Figure 3.5 shows <j> in a search region of f T , (t 0 ) = ±1° for three sample central 
transfer angles. Each crossing of n7t indicates a choice of maneuver point which is robust 
to bum magnitude error. Figure 3.6a shows the set of all such points for transfer angles 
in the range of 45° to 180°. There are two robust maneuver points for each transfer angle. 

When </) = ~ n, the condition for burn direction error tolerance (Equation 3.6) is 

satisfied. Unfortunately, there are no maneuver points that provide insensitivity to bum 
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Figure 3.6 a&b Bum Execution Error Tolerant Maneuver Points for Cocircular Staging 
Orbit 
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direction error (no — n crossings) for the circular orbit case studied. This means 
that there are no points at which a burn direction error is perfectly aligned with the 
allowable deviation in bum. There are regions, however, where these vectors are near to 
alignment. The component of the direction error along vv, t causes a slip in intercept time. 
The component which is perpendicular to vv (i causes a miss distance at intercept. If we 
minimize this perpendicular component (by choosing maneuver points where the vectors 
are nearly aligned), we will minimize the miss distance. 

For ignition time errors, we examine the angle 0 between the change in required bum 
caused by a delay in burn execution and the change in required bum corresponding to a 
delay in intercept. (Figure 3.7). 


ui,, 

Figure 3.7 Angle Between Change in Required Burn Caused by Execution Time Slip and 
That Corresponding to an Intercept Time Slip 

When <f> = rut, the condition for tolerance of ignition time errors (Equation 3.12) is 
satisfied. Figure 3.6b shows the set of all such points for transfer angles in the range of 
45* to 180*. There are two robust maneuver points for each initial chaser position 
throughout the range of transfer angles examined. 

Figure 3.8 shows the intercept time slips required to compensate for each type of bum 
execution error. As expected, an increase in applied burn magnitude will cause the chaser 
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Figure 3.8 Deviation in Transfer Time for Cocircular Staging Orbit 
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to arrive at intercept sooner (the vectors are anti-parallel: 0 is an odd multiple of 7t), and, 
in general, a delay in the bum will cause a delay in the intercept 


When using near-solution points to compensate for burn direction errors, the intercept 
time slip is 




SCAvyw,, 


(3.14) 


3.4 Application of Burn Error Tolerance Conditions to Highly Elliptic 
Orbits 

For highly elliptic orbits, we must choose f c (t 0 ), fy c (t 0 ), and 0 to specify a nominal 
rendezvous. 


t 

i 

i 



Figure 3.9 True Anomaly and Central Transfer Angle for Elliptic Orbits 

We will choose / c (f 0 ), and 0 and examine (p for a range of fy^o) (Figure 3.10). For 
elliptic orbits, each choice of / c (f 0 ) gives us a different solution, so we will fix 0 
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(0=1 35‘) and examine the entire range of initial chaser positions. 2 Figure 3.1 1 shows <p 
in a search region of fy (t 0 ) = ±5° for three initial chaser locations. 



Figure 3.10 Search Region for Elliptic Orbits 


Figure 3.12 a,b,&c show maneuver point conditions for transfer from the scaled elliptic 
staging orbit. Note that in addition to bum magnitude insensitivity points, direction error 
insensitivity points exist for this scaled elliptic case. Unlike the other regions, for the set 
of magnitude solutions that are enclosed by direction error solutions (near f c {t 0 ) = 120°), 
<P is an even multiple of 7t and an increase in the magnitude of the bum actually results in 
a delay in intercept. Figure 3.13 a,b,&c show the associated rendezvous time slips. The 
upper curve in Figure 3.12b causes extremely large time slip errors (Figure 3.13b), 
violating the linearity assumptions ( ^^ «1). These maneuver points should be 
avoided. Figures 3.14 and 3.15 do the same for transfer from the coelliptic staging orbit. 

Figure 3.16 describes the target-centered curvilinear coordinate system that will be used 
in the relative motion plots. 


2 For alternate solution sets, see Appendix C. 
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Figure 3.12 Burn Execution Error Tolerant Maneuver Points for Scaled Elliptic Staging 
Orbit and 135* Central Transfer Angle 
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Figure 3.13 Corresponding Deviation in Transfer Time - Scaled Elliptic Staging Orbit 
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Figure 3.14 Burn Execution Error Tolerant Maneuver Points for Coelliptic Staging Orbit 
and 135* Central Transfer Angle 
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Burn Magnitude Error 





Figure 3.15 Corresponding Deviation in Transfer Time — Coelliptic Staging Orbit 
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Figure 3.16 Curvilinear Coordinate System 


Figure 3.17 shows the relative motion plots tor transfer from a burn magnitude error 
insensitivity point (/ c (r 0 )= 45* and / % (f 0 )=+ 2.6379* from Figure 3.12a) for each of the 
three types of burn execution error. Note that in the case of the burn magnitude error, the 
three dispersed trajectories focus at the target (although at different times). The other two 
types of errors result in large miss distances. Figure 3. 1 8 repeats this for a bum direction 
error insensitivity point (/ c (r o )=90* and / % (t n )=-0.6416” from Figure 3.12b), and 
Figure 3.19 for a burn ignition time error insensitivity point (f c (t 0 )=90 and 
(jo) = +o 2738° from Figure 3.12c). Relatively large errors are used in all three cases 

so that the effects will be visible. 
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Figure 3.18 Relative Motion for Burn Direction Error Insensitivity Point 
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Chapter 4 

Navigation Error Tolerance 


Any error in the estimate of the state of the chaser, either position or velocity or a 
combination of the two, will result in an error in the estimate of the required bum. If the 
expected error in the estimate of position as well as the relationship between the expected 
error in the estimate of velocity and the expected error in position can be identified, the 
state error vector becomes one-dimensional. The freedom obtained by allowing a small 
variation in the intercept time can then be used to develop maneuver point conditions 
which are robust to such one-dimensional errors due to state uncertainties. 

4.1 Measurement Effects 

The errors in a vehicle's estimation of its state relative to the target are determined by the 
instruments making the measurements. In many cases the instrumentation may give 
direct information in only one direction, and information must be inferred in the other 
directions. Sometimes this leads to families of states which are indistinguishable from a 
navigation point of view, because they result in essentially the same measurement 
history. Figure 4.1a shows an example of states that are indistinguishable with angular 
measurements, and Figure 4.1b, states that are indistinguishable with range 
measurements, at least in a low dynamics environment. 
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a) Scaled Relative Geomedy b) Rotated Relative Geometry 


Figure 4.1 Indistinguishable Measurement Histories 

4.2 Scaled Relative Geometry Error 

During the early stages of rendezvous, the active vehicle may rely on "angles-only" 
navigation from star tracker measurements. The star tracker measurements supply direct 
information about the direction of the Line-of-Sight (LOS) from the chaser to the target. 
However, the chaser will have only very limited knowledge of its range along that LOS. 
Because of this, the relative state of the chaser will only be known within some scale 
factor. The error in estimated position of the chaser is primarily along the chaser’s 
relative position vector. The error in estimated velocity is primarily along the relative 
velocity vector and the ratio of the velocity error to the position error is the same as the 
ratio of the relative velocity to the relative position. As Figure 4.2 shows. 
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L*, = (! + £)£„, **8r = er ( 


(4.1) 


v«, = (l + £)y„, <*Sv = ev, 



Figure 4.2 Scaled Relative Geometry Error 


4.3 Rotated Relative Geometry Error 

The vehicle may, instead, rely on radar measurements. These measurements supply 
direct information about the range from the chaser to the target, and the chaser will have 
only very limited knowledge of the direction of the LOS to the target. For small 
rotations, the error in estimated position of the chaser is primarily perpendicular to the 
chaser’s relative position vector. Likewise, we may expect the error in estimated velocity 
to be perpendicular to the relative velocity. As Figure 4.3 shows, 

Lact = (/ + yx)r„, <=>5r = (ux)r„, 

Vac, = (! + <=> 5v = (tu)v„, (4.2) 
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Figure 4.3 Rotated Relative Geometry Error 


4.4 Covariance Analysis 

To demonstrate that the scaled relative geometry relationship in Equation 4. 1 does, in 
fact, occur, a covariance analysis was used to generate representative statistics of state 

Target 
t = t 0 



Figure 4.4 Position Error Ellipses After One Hour of Measurements 


56 



errors through one hour of optical measurements. The target is in the highly elliptic orbit 
around Mars and the chaser is in the scaled elliptic staging orbit The chaser starts with a 
large, diagonal covariance matrix (Figure 4.4). Correlations develop quickly, and Table 
4.1 lists the results at the end of one hour of optical measurements 1 . 



size and orientation of 
position error ellipse 
in r CIT frame 
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v~ /T fram 
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fc 
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b 

4 > 
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b 

♦ 

45 

45.9427 

1327 

9 

-0.01 

0.6840 

0.0051 

-0.82 

135 

135.4949 

2256 

26 

-0.02 

0.2154 

0.0288 

3.85 

180 

180.0474 

8848 

5 

-0.00 

0.3869 

- 

0.0167 

0.10 

225 

— ^ — 

225.0423 

— 

5003 

1 

0.00 

0..5134 

0..0023 

-0.02 

315 

315.0552 

1654 

1 

0.00 

0.8142 

0.0004 

0.00 


Table 4.1 Orientation of Position and Velocity Error Ellipses after One Hour Of Optical 
Navigation 


1 a is the semi-major axis, b is Ihe semi-minor axis, and <J> is the angle between the major axis and the 
chosen coordinate system . 
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As expected, the long major axis of the equi-probability position error ellipse is 
essentially parallel to the relative position vector and the much shorter minor axis is 
perpendicular to it. The long major axis of the velocity error ellipse is essentially parallel 
to the relative velocity vector and the much shorter minor axis is perpendicular to it. 
Now that the direction of the estimation error has been identified, a maneuver point 
condition robust to errors due to optical measurements can be derived. 

In contrast, the expected rotated relative geometry relationship of Equation 4.2 does not 
develop. The same covariance analysis was used to generate representative statistics of 
state errors through one hour of range measurements. Table 4.2 lists the results at the end 
of the hour of range measurements. 




size and orientation of 
position error ellipse 
in r c/T frame 

size and orientation of 
velocity error ellipse 
in yZ /r frame 

fc 
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0.0043 

-5.92 
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135.4949 

2277 

20 

-89.99 

0.4453 

0.0079 

-41.86 

180 

180.0474 

13739 

20 

90.00 

0.4317 

0.0456 

- 7.56 

225 

225.0423 

1302 

20 

89.99 

0..1147 

0..0167 

54.79 

315 

315.0552 

108 

19 

-87.36 

0.0704 

0.0076 

-53.56 


Table 4.2 Orientation of Position and Velocity Ellipses after One Hour Of Range 
Navigation 
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As expected, the position error is primarily perpendicular to the relative position vector . 2 
Strong correlations do develop in the components of the velocity error, but not the 
correlations expected of a rotated relative geometry. The orbital dynamics are creating 
more complicated indistinguishable states. Before robust maneuver point conditions can 
be developed for this case, more work needs to be done to understand the velocity error 
correlations. 

4.5 Scaled Relative Geometry Error - Burn Magnitude Error 
Approximation 

Error in the estimate of the position and velocity of the chaser will cause a post-bum 
velocity error which is, in general, a combination of bum magnitude error and bum 
direction error. Because we’re using staging orbits of nearly the same size and shape as 
the target orbits, the difference between the state of the chaser and that of the target is 
fairly small. If the chaser's state can be treated as a perturbation of the target's state. 


£c(*o) — £r(*o) 5r M 

v c (fo) = ¥r('o) + 5 v('o) 


the chaser's state at some later time can be written as 



= 

'8r(t 0 ) 


- k , *r,l 

'8r(t 0 ) 

_<5v(r,) 

K 

_8v(t 0 ) 


* 
1 

_Sv(t 0 )_ 


(4.4) 


^The cases which passed through regions of higher orbital dynamics ( f c - 45* and also f c -315) have 
gained significant information in both directions (a and b are both small). 
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K _ 

a: 2 =.M> 

dr{t 0 ) 

2 <?v(r 0 ) 

K _ <?v(f,) 

K - <?v(*i) 

3 dlih) 

4 Mt 0 ) 


In order for the chaser to intercept the target, the final deviation of the chaser position 
from the target position must be zero. Using this to solve for the deviation from target 
velocity required to achieve intercept, 

= -K^KfirUo) ( 4 . 5 ) 

The bum required for intercept is 

= Sv rtq - 8v(t 0 ) ( 4 . 6 ) 

In reality, the chaser’s relative position and velocity are only known within a scale factor. 
The actual position and velocity deviations are some small fraction larger or smaller than 
estimated: 

5r(/ 0 L, = (l + £)<5r(/ 0 )„, 

)ac, = 0 + £)Sv(t 0 K, ( 4 . 7 ) 

The actual required velocity deviation is 

=~K?KMt „)_ 

= -K;'R \(l + e)SrW„ 

= d+e)(^)„ (4.8) 
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and the actual bum required is 


( a ^L = ( 5 ^L~^ (r o ) - 

= (1 + e)(8v r ,') at - (1 + e)5v(t 0 )„, 

= (l + e)(Av„,) tw (4.9) 

Scaled relative geometry results in a scaled required velocity and a scaled required bum 
(Figure 4.5). The estimate of the required bum is correct in direction, but wrong in 
magnitude. This scaled relative geometry error can be treated as a bum magnitude error 
and the previously developed condition (Equation 3.3) applies. 



Figure 4.5 Scaled Relative Geometry Error Results in Bum Magnitude Error 


4.6 Scaled Relative Geometry Error — Maneuver Point Condition 


When the state of the chaser cannot be treated as a deviation from the state of the target, 
the previous applied method can be used to develop a more general form of the maneuver 
point condition which provides insensitivity to scaled relative geometry errors. 
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8 ( Av^ ) = X 0 SAt + L\ 8r c (t 0 ) + LiSr^)- 8vl (r 0 ) 


(2.3b) 


The error in the estimated position of the chaser is in the direction of the relative position 
vector with a magnitude which is some fraction of that of the relative position vector. 
The error in the estimated velocity is in the direction of the relative velocity vector and its 
magnitude the same fraction of the magnitude of the relative velocity. 

o) = £(£c(fo)-M'o)) = ££c/r(^o) 

8yr c (t 0 ) = e(vc (t 0 )-v T (t 0 )) = ev- c/T (t 0 ) (4. 1 0) 

Since transfer time and final position of the target do not change, the change in required 
bum due to a scaled relative geometry error is 

8(Av nq ) f(Z/ir c ^j-(r 0 ) (4.11a) 

For convenience, we define 

y¥.s = Leu (*o ) — Y.cn (*o ) 

8(Av nq ) = ew s (4.11b) 

In order for rendezvous to occur, this deviation in required bum combined with the 
deviation in required bum associated with a deviation in time of intercept must be zero. 


8(Av rtq ) = w h 5t x + £vv, = 0 

=-£w. 
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(4.12) 


The condition for insensitivity to a scaled relative geometry error is 


w, x vv (i = 0 


(4.13) 


Equation 4.13 is the maneuver point condition which provides insensitivity to scaled 


relative geometry error. 


The corresponding variation in time of intercept is 


8t x 


w. • w. 

— e — — = ±e 

w. • w . 

— 'l — *1 



(4.14) 


4.7 Deterministic Results 

Figure 4.6 shows the solutions for the bum magnitude error approximation and the more 
general scaled relative geometry maneuver point condition. In the case of the circular 
rendezvous, the results which are almost indistinguishable. However the two sets of 
results for the scaled elliptic and coelliptic rendezvous appear to be quite different 
(Figures 4.7 and 4.8). 

Even though the bum magnitude approximation is not exactly the same as the scaled 
relative geometry maneuver point condition, the difference is not necessarily significant. 
Because the scaled relative geometry error primarily causes a bum magnitude error, 
maneuver points which are insensitive to one type of error are fairly tolerant of the other 
type. 
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Figure 4.7 Comparison of Bum Magnitude Error and Scaled Relative Geometry Error 
Insensitivity Points for Scaled Elliptic Staging Orbit 
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Figure 4.8 Comparison of Bum Magnitude Error and Scaled Relative Geometry Error 
Insensitivity Points for Coelliptic Staging Orbit 
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The relative motion plots for transfer from one scaled relative geometry error insensitivity 
point ((/ c (f 0 )=90° and /^(f 0 )=1.07079° from Figure 4.7b) is shown in Figure 4.9. 

4.8 Monte Carlo Analysis 

The covariance matrices found in Section 4.4 were used to generate ten representative 
state dispersions around five navigation insensitivity points. These "actual" states were 
then propagated forward to rendezvous. Figures 4.11 through 4.15 show the relative 
motion plots from dispersed initial states for both optical and range measurements. 



Figure 4.10 Monte Carlo Analysis 

In the optical case, the pure scaled relative geometry portion of the error is fully 
compensated for by the appropriate choice of maneuver point. This leaves a much 
smaller initial bum error which causes a much smaller final position error. In Figure 
4.12, the dispersed states do not converge as well at rendezvous. This is the case in 
Section 4.4 where the scaled relative geometry had not fully developed during the one 
hour allotted. 
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Figure 4.9 Relative Motion Plots for Scaled Relative Geometry Insensitivity Point 


68 






y re , (km) Yrei (k m ) 




Figure 4.1 1 Dispersed States - / c (f 0 ) = 45* 
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Figure 4.13 Dispersed States - / c (r 0 )=180* 
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Figure 4.15 Dispersed States - / c (/ 0 )=315* 
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Chapter 5 
Conclusions 


5.1 Summary of Results 

The purpose of this thesis has been to develop maneuver point conditions which 
determine transfer phases robust to errors in burn execution and navigation. Chapter 2 
identified the degree of freedom gained by allowing a small variation in time of intercept 
and noted the restrictions placed upon use of this freedom by linearity assumptions. 
Chapters 3 and 4 showed that in cases where one dimensional errors in bum execution 
and navigational uncertainties can be aligned with the degree of freedom, the rendezvous 
will still occur, albeit at a slightly different time. Such insensitivity points can be applied 
to both circular and highly elliptic orbits and still maintain reasonable times of transfer 
and fuel consumption. In cases where the error of concern is never perfectly aligned with 
the degree of freedom, no maneuver points exist . However, in such cases, the maneuver 
points where the error and the allowable error come closest to alignment provide the most 
tolerance of that particular error. In general, for any one-dimensional error, the possible 
transfers can be searched for regions where the error and the allowable error come closest 
to alignment. These regions will provide a measure of robustness for that error. 
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5.2 Future Research 


This thesis developed 3 set of robust maneuver point conditions, and in doing so raised 
many questions. Can we establish conditions for the existence of solutions? If solutions 
do exist, how many exist? Why do we have three solution points in some regions for 
insensitivity to scaled relative geometry errors when there are only two or no solution 
points for insensitivity to any of the three bum execution errors? 

These maneuver points are only valid as long as the linearity assumptions are not 
violated. If the anticipated errors will necessitate time slips and vehicle motion beyond 
the region of linearity second order effects might be included to expand this region of 
linearity. 

This thesis applied the robust maneuver point conditions to circular and highly elliptic 
orbits. Further studies should investigate the range of orbit eccentricities and of central 
transfer angles. The relative size and orientation of the chaser orbit (angular difference in 
the line of apsides) will also effect the solution space. In addition, more complicated 
phasing orbits, such as those used in Shuttle missions should be investigated. The 
analysis will also vary slightly if the transition is to a station keeping point near, but away 
from, the target. This study could be expanded to direct rendezvous, i.e., no intermediate 
staging orbit. These conditions may even be applicable to rendezvous with a target in a 
hyperbolic orbit. 


Another area of research involves the identification of state uncertainties resulting from 
range measurements. While the position error is clearly perpendicular to the relative 
position vector, the relationship of the velocity error to this position error has not been 
identified. It may be that the dynamics of the orbit, coupled with the higher order terms, 
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may yield further information concerning this relationship. Any other expected state 
errors that develop while taking other types of measurements (or combinations of 
measurements) should also be investigated. 

Armed with this method, the mission planner can locate robust maneuver points that 
balance other mission constraints such as fuel expenditure, closing velocities, and final 
approach geometries. Out-of-plane errors should be examined in future studies. In 
addition, the perturbation effects due to central body oblateness, atmospheric drag, and 
the gravitational attraction of the Sun should not be ignored. Further research in these 
areas and others is required if the mission to Mars is to become a reality. 
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Appendix A 

Alternate Forms of the Maneuver Point Conditions 


A.1 Determination of Freedom in Final Relative Position 

The final position of the chaser is a function of the transfer time involved and the initial 
position and velocity of the chaser: 

r c (O = f(Ar,r c (r 0 ),v+(r 0 )) 


If there are deviations in any of these variables, the final position of the chaser changes. 


Sr At , ) = SAt + SrM + 

<?A t dr c (t 0 ) dv c (t 0 ) 


or more compactly. 


Sr c (t l ) = vl(t l )8At + K l 5r c (t 0 ) + K 2 dv + c {t ii ) 


Specifically, if we allow a small perturbation in total transfer time, 8t, the initial position 
and velocity of the active vehicle do not change. The final position deviations of the the 
active and passive vehicles are 
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££c(0 = Vc('i)&i 
Sr T (t 1 ) = v T (t l )8t l 


The difference between these two final position deviations is the deviation in final 
relative position of caused by the small slip in total transfer time 


Sr c (r, ) - Sr T (f, ) = (y^ (r, ) - y r (f, ))&, (A. 1 ) 

This is the direction of freedom in final relative position. Small errors in final relative 
position which are aligned with the relative velocity are absorbed as small variations in 
transfer time. Those which are perpendicular remain as pure miss distance. 

A.2 Magnitude and Direction Errors 

There is no deviation in initial position or pre-burn velocity of the chaser associated with 
a pure bum magnitude or direction error. 

An error in burn magnitude, <5(Ay)„ = eAv req , causes an error in the final position of the 
chaser: 


Sr c (0 = £K 2 Av rtq 


If this final position error happens to be aligned with the final relative velocity, the 
maneuver point provides some measure of insensitivity to this type of bum error: 

Kii&Yrtq') x {¥c(h )~ ir(^i )) = Q (A.2) 
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An in-plane error in bum direction, «(. Ay), s (Ay,„ x i h )S9, causes an error in the final 
position of the chaser: 

SccU,) = K 2 (Ay„,xi t )SS 

If dris happens to be aligned with the final relative velocity, the maneuver point provides 
some measure of insensitivity to this type of burn error. 

K 2 ( Ay„„ ) • (yc (f, ) - v r (r, )) = Q (A3) 

A.3 Ignition Time Error 

If we delay execution of the burn by some small amount of time, 8t 0 , the chaser s initial 
state, both position and pre-burn velocity, will change slightly. 

5r c (r 0 ) = yc(r 0 ) 5r o 

5yc(r 0 ) = £ c ('o)& 0 = -j 77 /- c(f o )&0 

Because of this delay, we have less time to arrive at the nominal intercept point: 

SAt = -8t 0 

Therefore, the deviation in the final position of the chaser is 
^Cc( r i) = [ _ Hc( r i) + ^1 ¥c( f o) + ^2lc^°^ f " 
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In order for rendezvous to occur, this deviation in final position combined with the 

deviation in final relative position associated with a deviation in time of intercept must be 
zero. 


[ Yc ('i ) + K \ Yc ('o ) + K ig c )]& 0 + [v“ ('i ) ~ Yr )]&, = 0 

For this to be true, the two vectors must be parallel or anti-parallel: 

[~Yc (', ) + Yc (t 0 ) + K 2 g c (t 0 ) ] x (v- (/, ) - v r (r, )) = 0 ( A 4) 

Equations A.2, A.3, and A.4 are the maneuver point conditions which provide 
insensitivity to burn magnitude, direction, and ignition time errors. They are equivalent 
to those developed in Chapter 3. 
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Appendix B 

Other Mission Constraints 


Once the valid maneuver points have been identified other mission constraints can be 
used to make the final choice of transfer. Fuel expenditure will always be an important 
constraint. Closing velocity should be minimized for safety and control. Of course, 
when fuel expenditure is constrained, closing velocity tends to follow. The final 
approach angle should be chosen so that the target is against a dark background above the 
chaser. Plots for these three constraints are presented below, but similar plots should be 
made for each constraint of importance and consulted when making the final choice of 
transfer conditions. 


See figures B.la-c for fuel expenditure, closing velocity, and approach angle 
corresponding to the insensitive maneuver points for the cocircular Earth orbits. 

See figures B.2a-c and B.3a-c for fuel expenditure, closing velocity, and approach angle 
corresponding to the coelliptic and the scaled elliptic staging orbits. 



Figure B.l LVLH Elevation Angle of LOS to Target 
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Figure B.2a Fuel Expenditure - Cocircular Staging Orbit 


84 








Av tot (m/sec) Av T0T (m/sec) Av T0T (m/sec) 







Burn D 

irection Error 


























120 180 240 300 360 



fc(to) (deg) 


Figure B.3a Fuel Expenditure - Scaled Elliptic Staging Orbit 
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Figure B.3c Final Elevation Angle of LOS to Target - Scaled Elliptic Staging Orbit 
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Figure B.4a Fuel Expenditure - Coelliptic Staging Orbit 
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Figure B.4b Closing Velocity - Coelliptic Staging Orbit 
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Appendix C 
Alternate Solution Sets 


This appendix contains two sets of scaled elliptic maneuver points for comparison. 

C.l Alternate Fixed Central Transfer Angle 

The first shows the maneuver points providing burn error tolerance for a 90 central 
transfer angle. Because less time is involved (than for 0 = 135”) these points are grouped 
much more closely. For larger transfer angles, the region where solutions exist will be 
larger. A fixed central angle results in different transfer times in different regions of the 

orbit. 

C.2 Fixed Transfer Time 

The second set shows the maneuver points providing burn error tolerance for a fixed 
transfer time of 3/8 of the target's orbital period. Transfers near periapse involve a larger 
central transfer angle, and those passing through apoapse, a smaller transfer angle. When 
planning rendezvous with a fixed transfer time, the resulting central transfer angle should 

be monitored. 
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Figure C.1 Maneuver Poims Providing Burn Execu.ion Error Tolerance for Scaled 
Elliptic Staging Orbit and 90‘ Central Transfer Angle 
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Figure C.2 Maneuver Points Providing Burn Execution Error Tolerance for Scaled 
Elliptic Staging Orbit and Transfer Time = 3/8 Pj 
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